Coherence and dynamics of a two-component condensate coupling with a continuum 
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We investigate the many-body coherence and population dynamics of a two-component atomic 
Bose-Einstein condensate coupled with a continuum. We reveal that the combination of dissipation 
and nonlinearity may induce different steady states and modify the bifurcation of coherence. The 
steady states depend on the relative dissipation strengths. The coherence may be even enhanced 
by the dissipation effects. The combination of dissipation and nonlinearity allows one to control the 
switching between different self-trapped states or the switching between a self-trapped state and a 
non-self-trapped state. This opens an alternative route for manipulating coherence and population 
transition in an open many-body system via dissipation and nonlinearity. 

PACS numbers; 03.75.Kk, 03.75.Lm, 05.30.Jp, 03.75.Gg 



I. INTRODUCTION 

Since the experimental realization of the two- 
component condensates, such as Bose-Einstein conden- 
sates (BECs) of in two different hyperfine levels [ll- 
0] , the dynamical features of two-component condensates 
have been widely investigated in recent years. Within 
the mean-field description, several macroscopic quantum 
behaviors and interesting dynamical properties 0413 
have been explored. Within the quantum-field descrip- 
tion, several many-body quantum phenomena [2ll429{ . 
such as quantum self-trapping ^21j and spontaneous sym- 
metry breaking [l^, have been predicted. Some of these 
dynamical properties have also been observed in experi- 
ments [il, . 

In real experiments, due to the atomic loss induced 
by inelastic collisions and the coexistence of condensed 
and non-condensed atoms [3ll - l33| . the atomic conden- 
sates are not absolutely closed systems. The dissipative 
processes can be effectively described by non-Hermitian 
(NH) Hamiltonians j34i-39]. The dissipation- induced ef- 
fects in a NH Bose-Hubbard dimer with a complex on- 
site energy (40l - |43| or a complex coupling term (4^ have 
been studied extensively. If a system has an intrinsic 
mechanism balancing the dissipations, the dissipation 
can lead to a constructive effect, such as the enhanced 
self-trapping [4l|, [4^ and the inhibited losses of atoms 
[45| . Moreover, the dissipation-induced coherence in an 
open two-mode BECs system have been studied by the 
master equation method [i^ H^l . 

In single-particle systems, the interaction between a 
two-state model and a continuum state is a paradigm 
to understand the crossing of quasi-stationary states and 
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its transition [48|, |49|. Up to now, there is still no any 
study on the coupling between a many-body two-state 
system with continuum. The experimental realization of 
two-component BECs provides a new possibility for ex- 
ploring many-body quantum phenomena of atoms cou- 
pled with continuum, such as many-body coherence, dis- 
sipative dynamics and population transition. In particu- 
lar, the study of the interaction between two-component 
condensates and continuum will provide practical appli- 
cations in manipulating the two-component condensates. 

In this article, we investigate the many-body quantum 
coherence and population dynamics in a two-component 
condensate coupled with a continuum. Our system is 
described by a NH Bose-Hubbard dimer whose diago- 
nal and off-diagonal elements are both complex numbers. 
In the mean-field theory, the system obeys a two-mode 
NH Gross-Pitaevskii Hamiltonian, in which the nonlin- 
ear terms describe the atom-atom interactions. We find 
that the combination of dissipation and nonlinearity may 
induce different steady states and modify the bifurca- 
tion of coherence. By tuning the nonlinearity and dis- 
sipation strength, it is possible to observe the coherence 
enhancement. Particularly, different relative dissipation 
strengths between two hyperfine levels will drive the two- 
component condensates into different stable states. In 
the time evolution, the system can jump from one sta- 
ble state to another stable state and the switching time 
depends on the nonlinearity strength. Our results show 
that the combination of dissipation and nonlinearity can 
be used to manipulate the steady behaviors, such as, con- 
trolling the transition between two self-trapping states 
or between a self-trapped state and a non-self-trapped 
state. Therefore, our results provide an alternative route 
for manipulating the many-body coherence and popula- 
tion dynamics by utilizing dissipation and nonlinearity. 

The structure of this article is as following. In sec- 
tion II, we give the non-Hermitian Hamiltonian and the 
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equations of motion. In section III, we analyze that how 
dissipation and nonhnearity affect the bifurcation of co- 
herence. In section IV, we show various nonhnear dissi- 
pative dynamics under different conditions. In the last 
section, we briefly summarize our results. 



II. NON-HERMITIAN TWO-MODE 
BOSE-HUBBARD MODEL 

We consider the interaction between a two-component 
condensate, whose atoms occupy two different hyperfine 
levels, and a continuum. If the trap potential is suffi- 
ciently strong, one can assume the condensed atoms only 
occupy the lowest vibrational state and so that the sys- 
tem obeys the Hamiltonian 



H = 



^ - (b\b2 + 4^1) + Ui2b\blb2bi 



E 

3=1.2 



b]b3 
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(1) 



Here &j (bj) are the bosonic creation (annihilation) op- 
erators for atoms occupying the j-th hyperfine level. In 
the complex coupling term, — i\/TiT2/'2, ft denotes 
the direct coupling between the two hyperfine levels and 
iy/TiT2/2 describes the effective coupling between them 
via the continuum [i^ . The inter-component inter- 
action Ujj and the intra-component interaction U12 can 
be controlled by adjusting the atomic s-wave scattering 
lengthes via Feshbach resonance techniques. In the com- 
plex energies, ej — iTj/'i, the imaginary parts Tj denote 
dissipation rates for atoms in the j-th hyperfine level due 
to the interaction with the continuum [3, El] . For the 
case of r2 — 0, the system (I) is just the model studied 
in [4II, 14^. For the case of no atom-atom interaction, 
Ujj = U12 = 0, the system (1) is reduced to the single- 
atom system [4^, . To simplify, we set fi, = I and adopt 
i7 as a unit to rescale the other parameters Ej, Tj, U12 
and Ujj. 

To investigate the dynamics of the system (I) in the 
Bloch representation, we introduce the angular momen- 
tum operators. 



L. = lib\b2 + b^bl), 

Ly = ^(fol62-fol4), 

= \iblb2-b\b^), 



(2) 



whose Casimir invariant is = {N/2){N/2 + 1). Here 
N = b\bi+b\b2 is the total number operator. With these 
angular momentum operators, the Hamiltonian (I) can 
be written as 



H ^2{n- iA)L^, + GLl + {S~ ij)L^ 



iTN, 



(3) 



where A = V^^/2, G = (C/n + U22 - 2U12), <5 = £2 - 
ei + iUii-U22)N/2, 7 = (r2--ri)/2 and T = (ri+r2)/4. 



To investigate the mean-field dynamics, introducing 
the SU{2) coherent states 

\xi,X2) = ^{xib\+X2bl)''\0,0), (4) 

with two complex coefficients Xj, the mean-field Bloch 
vectors Sfc — (L^.) /N, k — x,y, z, read as 

^ _ {xi,X2\Lx\xi,X2) _ {xlX2 + X1X2) 

{xi,X2\xi,X2)N 2n 

_ {xi,X2\Ly\xi,X2) _ {XIX2-X1X2) 



^ {xi,X2\xi,X2)N 2in 

{xi,X2\L^\xi,X2) _ {x'^X2-xlxi) 



{xi,X2\xi,X2)N 



2n 



(5) 



Here, n = + P is the norm. By using the method 
developed in Refs. |4ll, [l^j, the equations of motion for 
the Bloch vector are given as 



-iS + 2Cs,)sy 
(S + 2Cs^)s^ - 2ns 
2ns. 



y -r 4AS2SJ, - 7(1 



2js^Sa:- A{l-4sl), 

+ 4:ASxSy + 2jSzSy, 

~4sl)/2. (6) 



To obtain these equation, we have taken the semiclassical 
limit N ^ 00 with the mean-field interaction strength 
C — NG kept unchanged. Therefore, s"^ — s-^ + Sy + 
si = 1/4 is a constant and the dynamics are regular 
and confined onto the Bloch sphere. However, the total 
probability n decays as 



n = -2(2As^ + 75^ + T) 



(7) 



The Bloch equation (6) is equivalent to a generalized 
NH nonlinear Schrodinger equation. Introducing the un- 
normalized complex numbers ipj, which are associated 
with the coefficients Xj of the many-particle coherent 
state [4II, the Bloch equation (6) corresponds to the NH 
nonlinear Schrodinger equation 



'dt \il^2 



-S-CK-i{T -j) n^iA 

Q-iA 5 + CK-i{j + T) 



V'2 



(8) 



with K = (|'02p — + l^2p)- In recent years, 

similar NH nonlinear Schrodinger equations have been 
proposed to describe open atomic BECs [41 43, 50] and 
two-channel waveguides with gain and loss j51l | . 



III. COHERENCE BIFURCATION UNDER 
DISSIPATION AND NONLINEARITY 

In this section, we investigate the combined effects of 
the dissipation and nonhnearity on the phase coherence 
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between two modes. Here we use the contrast in interfer- 
ence experiments (52l . [53j to measure the phase coherence. 
The contrast is defined as 



2 I {b\b2 

{blh 



(9) 



In this work, we mainly focus on the coherence of the 
steady states (the fixed points of the equations of mo- 
tion). These steady states can be calculated numerically 
via Newton flow method from s = . For the case of 
(5 = 0, we have 

_ s°[Q(2ACs°-7f})-h4fl(A2-hl]2)] 



2[QAn + 2CsO(4A2(sO)2 + fiS)] 



Q 



(10) 



where Q = 7[1 — 4(s°)2] and denote the corresponding 
fixed points. By using Eq. (jlOl) and the normalization 
condition = (5°)^ -|- (§0)2 -|- (§0^2 _ ^^^^ determine 

the values of s° and then s° and s° from Eq. (fTO)) . A key 
issue is the stability of the fixed points. To do this, we 
linearize Eq. ^ around the fixed points by substituting 



Asx,Sy ^ s° + A.t 



As 



(11) 



with a small deviation (As^,, Asy, Asz) into Eq. ([6]), 
thereby one can obtain a linearized equation for this small 
deviation [H, IH, HF 



(12) 
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FIG. 1; (Color online) Dependence of the steady contrast a on 
the relative dissipation strength 7 and the imaginary coupling 
strength A. The nonlinearity C is chosen as (a) and (b) 1. 
(c) and (d) are sections of (a) and (b) for 7 = 1 and A = 1, 
respectively. The red solid (dotted) lines correspond to the 
stable (unstable) steady states for the nonlinear case, and the 
blue solid lines are steady states for the linear case. The other 
parameters are given as = 1 and 5 = 0. 



the solid lines (see Figs. 1(c) and 1(d)). In particular, we 
observe that the steady contrast a for the nonlinear case 
is larger than the one for the linear case. This means 
that the coherence can be enhanced by the combination 
of the dissipation and the nonlinearity. In addition, we 
find that the coherence may be enhanced by the imagi- 
nary coupling strength A, see Fig. 1(c). As A increases, 
the coherence slowly decreases first and then increases. 
This is similar to the dissipation induced coherence in 
previous works [i^ 113] . 



with the coefficient matrix 

275^ -I- SAfiO -2Cs53 
M= I 2Cs°-l-4AsO 4AsO-f27sO 
4As° 2n 



2^sl - 2C4 

2Csl -2n + 27s: 

4AsO + 475° 



(13) 



The fixed points are linearly stable if and only if there is 
no eigenvalues with a positive real part for the coefficient 
matrix M. The eigenvalues of the coefiicient matrix can 
be obtained numerically. 

We show our numerical results for the steady contrast 
a in Fig. 1. The dependence of a on the relative dissi- 
pation strength 7 and the imaginary coupling strength 
A are shown for two different cases: (a) the linear case 
C = and (b) the nonlinear case C = 1. In Figs. 1(c) 
and 1(d), we plot the sections of Figs. 1(a) and 1(b) for 
7 = 1 and A = 1, respectively. The other parameters 
are given as $7 = 1 and d — 0. The numerical results 
indicate that there exist two fixed points for the nonlin- 
ear case of C = 1 and while there only exists one fixed 
point for the linear case of C = 0. By implementing 
the linear stability analysis, we find that only one of two 
fixed points for the nonlinear case is stable, marked by 



FIG. 2: (Color online) The steady contrast a versus the non- 
linearity C for n = 1, 7 = 0.4 and A = 0.3. Here the red 
solid lines correspond to the stable fixed points, and the blue 
dotted lines correspond to the unstable fixed points. 

In Fig. 2, we show how the nonlinearity C affects the 
contrast a. For weak nonlinearity, the contrast a has 
two different values corresponding to one stable and one 
unstable fixed points. As the nonlinearity increases, a 
bifurcation appears at a certain critical value of the non- 
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FIG. 3: (Color online) Mean- field dynamics on the Bloch 
sphere for the Hermitian (A — 0, top) and non-Hermitian 
(A = 0.4, bottom) cases for Q = 1, S = and different values 
of the nonlinearity [C = (left), 2 (right)]. Different tincto- 
rial orbits on the Bloch sphere correspond to different initial 
conditions. 



linearity. This bifurcation of the contrast a is associated 
with the occurrence of one stable and one unstable points 
beyond the critical nonlinearity. Here the red solid lines 
denote the stable fixed points and the blue dotted lines 
represent the unstable fixed points. 



IV. NONLINEAR DISSIPATIVE DYNAMICS 

Based upon our understanding of the steady states in 
above, in this section, we analyze the mean-field dynami- 
cal behavior arising from the interplay of dissipation and 
nonlinearity. For the case that the dissipation strength 
in one of hyperfine levels is zero, i.e. A = 0, the cor- 
responding NH Hamiltonian is reduced to the ones in 
Refs. [411, . The mean- field dynamics has been studied 
in detail for this case, in which the dissi pat ion enhanced 
self-trapping states have been revealed ]4ll . [i^ . In this 
work, for simplicity, we focus on considering the following 
two cases: (A) 7 = and (B) 7 7^ and A 7^ 0. 





FIG. 4; (Color online) Mean-field dynamics on the Bloch sphere for the case 7 > A (top) and 7 < A (bottom) for 5 = and 
fl = l. (a)-(c) 7 = 0.4, A = 0.15; C = 0.0 (a), 1.5 (b), 5.0 (c). (d)-(f) 7 = 0.05, A = 0.324; C = 0.0 (d), 1.5 (e), 5.0 (/). Here 
different tinctorial orbits in the Bloch sphere correspond to different initial conditions with s(0) = [0, 0, 0.5] (red), [0.48, 0, 0.14] 
(blue), [0.3,0,0.4] (purple). 
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7 = 



7/0 and A / 



Under the condition of Fi = the relative dissipation 
between the two hyperfine levels vanishes, that is 7 = 0. 
Therefore, the mean-field Bloch equations are given as 



s'y ~ {S + 2Csz)sx — 20,3^ + 2AsxSy, 
with the total probability n satisfying 



A, 



h = - — (2.Sj, + f)?i. 



(14) 



(15) 



In Fig. 3, we show the dynamical behavior of the Bloch 
equation (fT4)) on the Bloch sphere for the Hermitian case 
of A = (top) and the non-Hermitian case of A = 0.4 
(bottom) with il = 1 and (5 = 0. For the Hermitian 
case of A = 0, the Bloch vectors evolve periodically on 
the surface of the Bloch sphere and form closed orbits 
dependent upon initial conditions. There are two stable 
centers at Sj, = = and Sx — ±1/2, which are shown 
in Fig. 3(a). When the nonlinearity strength C increases, 
one of the two centers becomes unstable and a saddle and 
two stable centers appear after a bifurcation, which are 
shown in Fig. 3(b). In the case of A = 0.4 (bottom), 
we see a drastic modification of these patterns. On the 
one hand, all orbits on the surface of the Bloch sphere be- 
come non-closed. On the other hand, the system starting 
from different initial states may relax to the same state 
of Sy = §2 = and = —1/2, see Figs. 3(c) and (d). 
Furthermore, the change of nonlinearity and dissipation 
only modifies the evolution paths and dissipation veloc- 
ities to the final state. Therefore, in the case of 7 = 0, 
for different nonlinearities and initial atom populations, 
the final state will be completely definite, that is, it is an 
equal population state. 

To better understand the mean-field dynamics, we an- 
alyze the fixed points and their stability. The fixed points 
are determined by 



(4) 



f72 -2Cf7sO -t-A2(sO)2' 

(sg)2j2[^2^A2(sg)2] 

+ [l]2_2OT4+A2(,sS)^)] 



1 

4 



0. (16) 



For the case of 5 



and 



0, we have two stationary states 
— ■5° = from Eq. However, 
considering Eq. ([T5|) together, we only have one stable 
fixed point 



(17) 



This explains why all initial states finally decay into this 
stable fixed point. 



Now, we consider the general case of 7 7^ and A ^ 0, 
which obeys the Bloch equations Eq. (6). In Fig. 4, we 
show the mean-field dynamics on the Bloch sphere for 
the two different cases: 7 > A (top) and 7 < A (bottom) 
with fl = 1 and (5 = 0. Similarly, the orbits on the surface 
of the Bloch sphere are not closed. In the case of 7 > A, 
it is observed that, for linear and weakly nonlinear cases, 
all different initial states will always evolve into the same 
final state, as shown in Figs. 4(a) and (b). However, for 
strongly nonlinear cases, different initial states will evolve 
into different final states. As shown in Fig. 4(c), for the 
case of C = 5, two different initial states evolve into dif- 
ferent final states. Therefore, for 7 > A, the change of 
the nonlinearity and initial states not only can modify 
the evolution paths and dissipation velocities to the final 
state but also may change the final state. However, for 
7 < A, independent upon the initial state and nonlinear- 
ity strength, the system will always evolve to the same 
point on the Bloch sphere, which corresponds to the same 
final state, while the change of nonlinearity modifies its 
evolution path and dissipation velocity to the final state, 
as shown in Figs. 4(d)-(f). After a further calculation, 
we find that for 7 > A, the system has two stable fixed 
points, while 7 < A, the system has only one stable fixed 
points. So the competition between 7 and A results in 
different final state. It is important to note that the dis- 
sipation can be controlled by shining a laser beam onto 
the condensates [56| and relative dissipation rate 7 can 
also be changed at the same time. Therefore the com- 
bination of nonlinearity and dissipation can be used for 
controlling the dynamics. 




5 10 15 20 t 25 30 35 40 



FIG. 5: (Color online) Time evolution of Sz{t) in the case 
of 7 = 0. Here (5 = 0, = 1, A = 0.2 and the initial state 
locates at the north pole. The blue and red curves correspond 
to C = 2.1 and C = 2.5, respectively. 

Below, we discuss how to control the population 
switching by controlling the nonlinearity and dissipation. 
We assume that all atoms initially occupy the second hy- 
perfine level, which can be easily prepared in labs [52| . 
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FIG. 6: (Color online) Time evolution of s^{t) in the case of 
7 > A. Here (5 = 0, = 1, 7 = 0.4, A = 0.15 and the initial 
state locates at the north pole. The blue, black and red curves 
correspond to C = 1.5, C = 2 and C = 5, respectively. 
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FIG. 7: (Color online) Time evolution of Sz{t) in the case 
of 7 < A. Here 5 = 0, O = 1, 7 = 0.05, A = 0.324 and 
the initial state locates at the north pole. The blue and red 
curves correspond to C = 1.5 and C = 5, respectively. The 
inset gives the enlarged region between the two dotted lines. 

This initial state corresponds to the north pole of the 
Bloch sphere. In Fig. 5, we show the time evolution of 
Sz (t) for the case of 7 = with two different values of the 
nonhnearity. In this case, it is found that an initial self- 
trapped state can finally evolve into an equal-population 
state. In Fig. 6, we show the numerical results for Sz(t) 
in the case of 7 > A with three different values of nonhn- 
earity. From this numerical calculations, we find that the 
system may switch from an self-trapped state to another 
self-trapped state. The numerical data show that strong 
nonhnearity C > 2 drives the system close to self-trapped 
state with large population imbalance, Sz ~ —0.5. In ad- 
dition, for 7 < A, the system will jump from an initial 
self-trapped state to an equal-population state, as illus- 



trated in Fig. 7. By varying the relative dissipation 
rate 7, one can observe the switching between different 
self-trapped states (see Fig. 6) or the switching from a 
self-trapped state to an equal-population state (see Figs. 
5 and 7). On the contrary, one can also understand the 
dissipative mechanism of an open two-component BECs 
system by measuring the population difference between 
two hyperfine levels [52i] . 

V. CONCLUSIONS 

In summary, we explore the many-body quantum co- 
herence and population dynamics in a two-component 
condensate coupled with a continuum. The system is 
described by a NH Bose-Hubbard dimer whose diagonal 
and off-diagonal Hamiltonian elements are both complex 
numbers. It is revealed that the combination of dissipa- 
tion and nonhnearity may induce different steady states 
and modify the bifurcation of coherence. By tuning 
the nonhnearity and dissipation, the coherence enhance- 
ment exhibits. Particularly, different relative dissipation 
strengths between two hyperfine levels 7 will drive the 
two-component condensate into different stable states. 
Under the condition of 7 = 0, the atoms always evolve 
into a balanced state with equal population. In the case 
of 7 7^ 0, dependent upon the values of 7 and A, the 
system will evolve into different final states. For 7 > A, 
the system always evolve into a steady state with self- 
trapping. The change of nonhnearity and initial state 
not only can modify the evolution path and the dissipa- 
tion velocity to the final state but also can change the 
strength of self-trapping. For 7 < A, the atoms always 
evolve to the same quasi-equal-population state with a 
small Sz in order of 10~^, although the evolution path 
might be different for different nonhnearity strengthes. 

Our results show that the combination of dissipation 
and nonhnearity can be used to manipulate the steady 
behaviors, such as, controlling the transition between two 
self-trapped states or between a self-trapped state and 
a non-sclf-trapped state. Therefore, our results provide 
an alternative route for manipulating the many-body co- 
herence and population dynamics by utilizing dissipation 
and nonhnearity. Based on the experimental techniques 
in Ref. it is possible to observe the results in this 
article. 
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